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Abstract 

We consider Brans-Dicke theory with a self-interacting potential in Einstein conformal 
frame. We introduce a class of solutions in which an accelerating expansion is possible 
in a spatially flat universe for positive and large values of the Brans-Dicke parameter 
consistent with local gravity experiments. In this Einstein frame formulation, the theory 
appears as an interacting quintessence model in which the interaction term is given by the 
conformal transformation. In such an interacting model, we shall show that the solutions 
lead simultaneously to a constant ratio of energy densities of matter and the scalar field. 
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1 Introduction 



Cosmological observations on expansion history of the universe indicate that the universe is 
in a phase of accelerated expansion (Riess et al. 1998 ; Perlmutter et al. 1997a,b). This 
phenomenon may be interpreted as evidence either for existence of some exotic matter com- 
ponents or for modification of the gravitational theory. In the first route of interpretation one 
can take a mysterious cosmic fluid with sufficiently large and negative pressure, dubbed dark 
energy. This viewpoint also includes quintessence models (Ratra and Peebles 1988; Caldwell, 
Dave and Steinhardt 1998; Steinhardt, Wang and Zlatev 1999; Macorra and Picinelli 2000; 
Gonzalez-Diaz 2000; Bludman and Roos 2002; Cai, Saridakis, Setare and Xia 2010), which 
invokes a scalar field minimally coupled to gravity. The scalar field takes negative pressure 
during its evolution by rolling down a proper potential. In the second route, however, one 
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attributes the accelerating expansion to a modification of general relativity. This includes 
scalar-tensor theories, scalar fields non-minimally coupled to gravity (Faraoni f999; Faraoni 

2004) . 

Here we shall consider a self-interacting Brans-Dicke (BD) theory (Brans and Dicke 1961) 
as a prototype of scalar-tensor theories. The original motivation of the BD theory was the 
search for a theory containing Mach's principle which has found a limited expression in general 
relativity. As the simplest and best-studied generalization of general relativity it is natural 
to think about the BD scalar field as a possible candidate for producing cosmic acceleration 
without invoking a quintessence field or exotic matter systems (Setare 2007; Setare and Jamil 
2010). In fact, there have been many attempts to show that BD model can potentially explain 
the cosmic acceleration. It is shown that this theory can actually produce a non-decelerating 
expansion for low negative values of the BD parameter uu (Banerjee and Pavon 2001). This 
conflicts with the lower bound imposed on this parameter by solar system experiments (Will 
1993; Will 2005). Some authors propose modifications of the BD model such as introducing 
some potential functions for the scalar field (Bertolami and Martins 2000; Mak and Harko 

2005) , or considering a field-dependent BD parameter (Chakraborty and Debnath), without 
resolving this problem. 

All the works in this context use Jordan frame representation of BD theory. It is however 
well-known that this theory, like any other scalar-tensor theories, can be represented in the 
so-called Einstein frame by using a conformal transformation (Faraoni 1999; Faraoni 2004). 
Although these two conformal frames are mathematically equivalent there are some debates on 
their physical equivalence. Here there is a point which should be made clear. Despite the fact 
that the problem of physical status of the two conformal frames is open there is a tendency 
in the literature to ignore this problem and to work in Jordan conformal frame. The reason 
may be related to reluctance in accepting of the violation of weak equivalence principle due to 
anomalous coupling of the scalar field to matter systems in Einstein frame. It is however im- 
portant to keep in mind that the physical metric should be singled out already in the vacuum 
sector of the theory and the coupling of a given metric to matter systems is determined by the 
physical significance ascribed to it (Magnano and Sokolowski 1994). Thus a criterion based on 
the coupling of matter with gravity would be effective only if the physical frame were deter- 
mined on an independent ground. Apart from this point, anomalous gravitational coupling in 
Einstein conformal frame does not necessarily mean violation of weak equivalence principle. 
There is still a possibility that the effective mass of the scalar field be scale dependent. In this 
chameleon mechanism (Khoury and Weltman 2004 a,b), the scalar field may acquire a large 
effective mass in Solar System scale so that it hides local experiments while at cosmological 
scales it is effectively light and can provide an appropriate cosmological behavior. 

Along these lines, we would like to consider Einstein frame formulation of the theory as a 
representation which provides different possibilities in a cosmological setting with respect to 
the Jordan frame. In particular, we deal with two important issues. In section 2, we focus 
on the question that whether it is possible to achieve accelerating expansion of the universe 
for sufficiently large values of the parameter uj. We will show that the answer is affirmative 
within a class of solutions which corresponds to a specific form of the potential function of 
the BD scalar field. In section 3, we consider the coincidence problem which concerns with 
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the coincidence between the observed vacuum energy density and the current matter density. 
Recently, there is a large amount of interest to realize this problem as a consequence of an 
interaction between matter systems and the dark sector (Zimdahl, Pavon and Chimento 2001; 
Zimdahl and Pavon 2003; Campo, Herrera and Pavon 2009; Bisabr 2010a). Although the 
whole idea seems to be promising, however, the suggested interaction terms are usually phe- 
nomenological and are not generated by a fundamental theory. In the model presented here, 
there is an interaction between the BD scalar field and matter systems induced by conformal 
transformations. We will investigate the consequences of this interaction term and derive an 
expression for the relative energy densities. We also study the evolution of the latter and show 
that even though the two types of energy densities evolve differently as the universe expands, 
the ratio can take a small constant configuration at late times. 



2 Accelerating Expansion 

We begin with a modified form of the BD action in Jordan frame 

S JF = J d 4 x^==g{4>R - ^V^V.0 - V (ft) + S m (g^, (1) 

where is the self- interacting BD scalar field with a potential function V(<f>), u is a constant 
parameter and S m is the action of matter which depends on the metric g^ v and some matter 
fields collectively denoted by ip. A conformal transformation 

9v* ->• fiV = ^ 2 9^u (2) 

with Q = y/G(p brings the above action into the Einstein frame (Faraoni 1999; Faraoni 2004). 
Then a scalar field redefinition 



with 0o ~ G^ 1 , (j) > and u > — § transforms the kinetic term of the scalar field into a 
canonical form. In terms of the variables (g^, y?) the BD action in the Einstein frame is 
(Faraoni 1999; Faraoni 2004) 

R 1 



S EF = J d 4 x^g-(— - -sTVwVvtp - U(<p)) + S m ( 9lw , 1>) (4) 



where 



Here V M is the covariant derivative of the rescaled metric g^ u and 



S m = J d A Xy/=g exp(-8J j y?) L m (g^^) (5) 



U(<p) = V(4>(<p)) exp(-8 J ^^V) (6) 
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is the Einstein frame potential. 

Variation of the action (4) with respect to g^ v and ip leads to the following field equations 



G, v = SnG(T, u + Tp (7) 

dUUp) 1 m 

Dip ±Zl = - a T 8 

Y dip 2 w 

where 

T*, = V^V^ - ^V 7 ^VV - C/(^)^ , (9) 



a = y and T = g^ v T^ u is the trace of the matter stress-tensor. Note that the parameter 
a is related to inverse of the BD parameter oj. In Einstein frame, the vacuum sector of 
the action consists of a scalar field minimally coupled to Einstein's gravity. The important 
difference between the Einstein frame representation of BD model and minimally coupled 
scalar field models is that in the former the scalar field interacts with matter systems. This 
anomalous gravitational coupling has no counterpart in Einstein's gravity. It implies that the 
stress-tensors of matter and the scalar field are not separately conserved. This can be easily 
checked by applying Bianchi's identities to (7) which leads to 

V%, = - V7£ = \a T V v ip (10) 

The parameter a measures the strength of the interaction. Here there are two important 
points in order : First, the parameter a is positive a > 0. It implies that energy transfer 
is from scalar field ip to matter systems. This feature is consistent with the second law of 
thermodynamics (Pavon and Wang 2009). Second, the model (4) should be constrained by 
local gravity experiments to avoid violation of weak equivalence principle. It is well-known 
that these constraints are satisfied when u » 1 (Will 1993; Will 2055). In the model (4) this 
translates into a « 1. This means that the theory can pass local tests if interaction of the 
scalar field p with matter fields is sufficiently small. We will return to this issue later. 
We apply the field equations (7) and (8) to a spatially flat Friedmann-Robertson- Walker space- 
time 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) (11) 

with a(t) being the scale factor. To do this, we take the matter system to be a pressureless 
perfect fluid (dust) with energy density p m . In this case, the gravitational equations (7) give 

d 2 

S- = k{p m + p v ) (12) 
cr 

r,d a 2 , 

2- + — = -kp v 13 
a a 2 

where k = 8irG, p^ = \ip 2 + U(ip) and = \ip 2 — U(ip). We may use the first equation to 
rewrite the second one as 

*+2- 2 = l -kp m + kU{p) (14) 
a cr 2 
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The equation (8) gives 



a dip 2 
On the other hand, the conservation equations (10) become 



(p + 3-v? + — f^- = --a p m (15) 



p m + 3-p m = Q (16) 



a 



p\ + 3-(w v + l)p v = -Q (17) 

(X 

where 

Q = -a (p p m (18) 

and uj^ = Pip/ p v is the equation of state parameter of the scalar field <p. The equation (16) 
can be solved which gives the following solution 

p m a 3 = Pmoe^ (19) 

where p m o is the present matter energy density in the universe. Now we introduce a class of 
solutions of the field equations, characterized by, 

V?=4= lna ( 20 ) 
v k 

in which j3 is a positive constant parameter of order of unity. If we put the latter into (15), 
we obtain 

- + 2- = --5- — ^ + ~®p m ) 21 
a p v? 2 

Comparing (21) with (14) leads to a consistency relation 

^ + J~ mM = -f(lf + (22) 

in which e = af3/2\fk = f3[2{2uj + 3)]~5 > 0. This consistency relation will be satisfied for 
an appropriate potential function U(ip). To find the form of this potential, we solve this first 
order differential equation which gives the following solution 

U(tp) = -l Pm ,ef^ + Ce-^o (23) 

-0 2 +e 

where 7 = /3 2 2+£ _ 3 and C is an integration constant. As a consequence, the relation (20) is 
a solution of the field equations for an exponential potential of the form (23). This double 
exponential potential is similar to the potential which is used in some quintessence models 
(Barreiro, Copeland and Nunes 2000; Sen and Sethi 2002; Pavlov, Rubano, Sazhin and Scud- 
ellaro 2002; Demianski, Piedipalumbo, Rubano and Tortora 2005). It is shown that this kind 
of potential of the quintessence field can lead to solutions which first enter a period of scaling 
through the radiation and matter domination eras and then smoothly evolve to dominate the 
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energy density for a wide range of initial conditions of the field"!" (Barreiro, Copeland and Nunes 
2000). Moreover, single exponential potentials are popular in modified f(R) gravity models 
(Sotiriou and Faraoni 2010; De Felice and Tsujikawa 2010). These models are conformally 
equivalent to BD models with potentials which their forms are closely related to the functional 
form of the f(R) functions (Bisabr 2009). In that context, single exponential potentials cor- 
respond to power law f(R) gravity models (Capozziello, Cardone, Piedipalumbo and Rubano 
2006; Bisabr 2010b). 

It is interesting to note that the relation (20) brings the law of evolution of matter density, 
given by (19), into the following form 

p m = p m0 a~ 3+£ (24) 

This is similar to the rule presented by some authors for characterizing decaying law of vacuum 
energy into dark matter (Wang and Meng 2005; Alcaniz and Lima 2005; Costa, Alcaniz and 
Maia 2008; Jesus, Santos, Alcaniz and Lima 2008; Costa and Alcaniz 2010). It states that the 
scalar field ip is constantly decaying into the matter so that the latter will dilute more slowly 
compared to its standard evolution p m oc a -3 . Since the observational lower bound imposed 
by solar system experiments on the BD parameter is u » 1, we should have e << 1 which 
means that evolution of matter density has a small deviation with respect to the standard one 
in Einstein's gravity. 

In the expression (23), the integration constant C can be determined by noting the fact that 
when we set = 0o ~ G^ 1 in the action (1), then V(<f>) characterizes the vacuum energy 
density corresponding to a cosmological constant, namely V{4>) = A/G. In this case, if — 
and then 

U(<p = 0) = V(<f>(ip)) = p^o (25) 
with p^o being the vacuum energy density in the Einstein frame. Applying the latter condition 
to the relation (23) gives C = p v0 + 7p m o- The potential function takes then the form 

U{<p) = l Pm o{e-^ - ef^) + p^e-^ (26) 
For this potential function, the Friedmann equation (12) becomes 

TH = t\t2 - ^moa- 3+£ + (^o + lV m0 )a-e 2 ] (27) 

where Q m0 = p m o/p c , Q<po — P^o/Pc and p c = 3iJg//c is the critical density. From the equations 
(13) and (27), it is straightforward to show that the deceleration parameter 



takes the following form 



1,,. 1 „Os 1 



*w - + -/') - (3 - - 2 m - 7 _ ( ^ +7)( , +1) _ 3+ . +g J -'} 

where we have used a(z) — (z + l) -1 . This relation gives deceleration parameter in terms of 
the redshift and constant parameters u and f3. We plot q(z) in fig. la for u = 40000. The 



1"Note that Einstein frame representation of BD models arc effectively equivalent to the so-called coupled 
quintessence models in which the quintessence field interacts with matter sector. 



6 




Figure 1: The plot of deceleration parameter q(z) (panel a) and the effective equation of 
state parameter (panel b), given by (29) and (31), for u = 40000 and some values of the 
parameter (3. We have set = |. The values j3 = 0.4,0.8,1 respectively correspond to 
£ = 0.0010,0.0020,0.0025. *° 



figure indicates that accelerating expansion is possible for positive and large BD parameter 
(u » 1) consistent with observations. Moreover, it shows that the signature flip of q corre- 
sponding to transition from decelerating to accelerating expansion takes place in recent times 
as experimentally demonstrated by (Riess 2001). 

We can also obtain the effective equation of state parameter, defined by, 

■"» = ^ (30) 

Combining the latter with (12), (13) and (28) leads to 

7e// = ^(2<7-l) (31) 

Evolution of j e ff with respect to the redshift is plotted in fig. lb. The figure shows that the 
BD scalar field does not appear as a phantom field in the Einstein frame since the curves do 
not cross the boundary 7 e jj = — 1. 



3 The Cosmic Coincidence 

One of the important features of the cosmological constant problem is the present coincidence 
between dark energy and dark matter energy densities (Steinhardt 1997; Zlatev, Wang and 
Steinhardt 1999). There is a class of models in which this observation is related to some kinds 
of interaction between the two components (Zimdahl, Pavon and Chimento 2001; Zimdahl 
and Pavon 2003; Campo, Herrera and Pavon 2009; Wetterich 1995; Amendola 2000; Tocchini- 
Valentini and Amendola 2002; Boehmer, Caldera-Cabral, Lazkoz and Maartens 2008). In these 
models the two components are not separately conserved and there is a flow of energy from 
dark energy to dark matter or vice versa. In this sense, dark energy and dark matter energy 
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densities may have the same scaling at late times due to the interaction, although they decrease 
with the expansion of the universe at different rates. The important task in this context is 
to find a constant ratio of dark energy to dark matter energy densities for an appropriate 
interaction term. Despite the fact that this approach seems to be promising, there is not still 
a compelling form of interaction which is introduced by a fundamental theory. Therefore one 
usually uses different interaction terms and tries to adapt them with recent observations. In 
BD theory presented in Einstein frame, there is a fixed interaction between the scalar field and 
matter sector induced by conformal transformation. The strategy that we are going to pursue 
in this section is to use this interaction term for finding evolution of the ratio r = p m / p^- 
To do this, we first write the expression (24) in the following form 

Pm = ^ m0 a" 3+£ (32) 

Then, we use (20) to write 

Pv = \f 2 + U(<p) = + STOnoM 2 - ^ 3+£ ) + 3%oa-? 2 } (33) 

We may combine the latter with (27) to obtain 

Pv = ^ {3 n m0 ^-^a-^ + ^^(70no + %o)a-? 2 } (34) 
Finally, expressions (32) and (34) give the ratio of energy densities 

r (z) = 3n m0 (z + i)^{3n m0 ^p(z + 1) 3 - £ + ^L( 7 ru + n^z + if}- 1 (35) 

We plot r(z) in fig. 2. The figure indicates that the ratio r(z) varies during expansion of the 
universe and tends then to some constant of order unity in recent times. This behavior seems 
to be independent of the value attributed to the parameter (5. However, the latter affect 
the variation of r at large redshifts. When (5 increases, which corresponds to increasing of 
e, variation of r(z) decreases. This implies that the effect of interaction of the two energy 
components is to make the ratio r have more smooth variations at large redshifts. 



4 Conclusion 

We have considered the possibility that a self-interacting Brans-Dicke field accounts for the 
accelerated expansion of the Universe in Einstein conformal frame. Our analysis is limited 
to a class of solutions characterized by (20) which is given in terms of a double exponential 
potential of the form (23). The first important feature of the solutions is that they provide a 
late-time accelerating expansion of the universe for large values of the BD parameter u. At the 
same time, these solutions modify the evolution of matter density to (24) which is the simplest 
possible way of stating that the matter dilution is attenuated due to its interaction with the 
scalar field ip. This evolution law indicates that the deviation from the standard evolution is 
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characterized by a positive constant parameter e which quantifies the decay rate. Accelerating 
expansion is possible for e « 1 (or u> » 1) in accord with local gravity tests. 
Another important feature of the model presented here is that it simultaneously solves the 
coincidence problem. As previously stated, the interaction of the BD field <p with dark sector 
has been recently taken as a natural guidance for addressing coincidence problem by some 
authors. However, in absence of an interaction or coupling term based on a fundamental 
theory, most of the current investigations have been limited to a phenomenological level. In our 
analysis, this interaction term is given by the conformal transformation. We have investigated 
the consequences of this interaction term and derived a relation giving the evolution of the 
ratio of energy densities r. We have found that in spite of the fact that the ratio can be 
variable during expansion of the universe in the past, it tends to some constant value of order 
unity at late-time thus solving the coincidence problem. 
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